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ABSTRACT 



Aims. To investigate the universality hypothesis of the initial mass function in the substellar regime, the population of the p Ophiuchi 
molecular cloud is analysed by including a new sample of low-mass spectroscopically confirmed members. To that end, we have 
conducted a large spectroscopic follow-up of young substellar candidates uncovered in our previous photometric survey. 
Methods. The spectral types and extinction were derived for a newly found population of substellar objects, and its masses estimated 
by comparison to evolutionary models. A thoroughly literature search was conducted to provide an up-to-date census of the cluster, 
which was then used to derive the luminosity and mass functions, as well as the ratio of brown dwarfs to stars in the cluster. These 
results were then compared to other young clusters. 

Results. It is shown that the study of the substellar population of the p Ophiuchi molecular cloud is hampered only by the high 
extinction in the cluster ruling out an apparent paucity of brown dwarfs. The discovery of 16 new members of p Ophiuchi, 13 of them 
in the substellar regime, reveals the low-mass end of its population and shows the success of our photometric candidate selection with 
the WIRCam survey. The study of the brown dwarf population of the cluster reveals a high disk fraction of 76^%. 
Conclusions. Taking the characteristic peak mass of the derived mass function and the ratio of brown dwarfs to stars into account, we 
conclude that the mass function of p Ophiuchi is similar to other nearby young clusters. 

Key words, stars: formation - stars: low-mass, brown-dwarf - stars: planetary system 
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1. Introduction 

Although there is a generally accepted theory of how a star 
forms through the collapse of a mo lecular cloud cor e and con- 
tracts to the main sequence (e.g., IShu et al ] l!987l) . two ma- 
jor breakthroughs in the mid-90s have added a new complex- 
ity to the star formati on field: the disco very of exoplanets 
dMavor & Ouelo dl 995l) and brown dwarf s (Qppenh eimer et all 
119951: iRebolo et al.lll995b iNakaiimaet aT] [l995)". Brown dwarfs 
are a key element in understanding the origin of stellar masses 
and their distribution, the initial mass function (IMF), and in 
identifying the plethora of hospitable hosts where planetary sys- 
tems can form. However, there is no current consensus regard- 
ing their formation mechanism and several scenarios are cur- 
rently proposed. Some defend a star-like formation, i.e. fragmen- 
tation of molecular clouds into low-mass cores driven by turbu- 
lence ( Pad oan & Nordlundll2002t iHennebelle & Chabrierl feOQS ) 
or gravity (Bonnell et al. 2008). But distinct formation processes 
are al so likely to take place, such as gr avitational instabilities in 
disks ( Stamatel los & Whitwo rth 2009), premature ejection from 
prestel lar cores (Re ipurth & Clarke! 1200 ll) . or photo-erosion of 
cores (iKroupa & Bouvien 120031) . The observational properties 
of brown dwarfs in young clusters show a global scaling-down 
trend from those o f stars, arguing in fav our of a common forma- 
tion scenario (e.g.. lLuhman et aO 2007). However, most studies 
in young star-forming regions suffer from incompleteness both 

* Based on observations made at the ESO La Silla and Paranal 
Observatory under program 085.C-0960. 



at lower masses (below 40 Mj up ) and in spatial content, fre- 
quently focusing on the inner regions of clusters. It is therefore 
unknown whether, as one moves to lower masses other forma- 
tion mechanisms take over (e.g., Whitwort h et al.ll2.Q10h . At the 
same time, the increasing capacity of numerical simulations to 
reproduce the collapse of entire molecular clouds is such that 
the comparison between simulations and observations starts to 
be hampered by our empirical knowledge based on small number 
statistics o f the low-mass end of the brown dwarf regime (e.g., 
lii[2009). Extensive studies of the lowest mass populations in 
different environments and stages of evolution are therefore re- 
quired to address these questions. 

To that end, several large photometric surveys are taking 
place to identify candidate brown dwarfs in various young 
cluste rs: for example, th e Galactic Cluster Survey (GCS) 
(e.g., lLodieu et alj 2011a) part of th e Infrared Deep Sky 



Survey (UKIDSS. lLawrence et al.ll2007l) at the UKIRT (United 
Kingdom Infrared Telescope), the SONY C survey (Substella r 
Objects in Nearby Young Clusters) (e.g. IScholz et al. 2009), 
or our WIRCam/CFHT ( Wide Field IR Cam era at the Canada 
France Hawaii Telescope, Puget et al. 2004) survey of nearby 
star-forming regions (e.g. iBurgess et alJ I200 9 J). among other 
substellar studies of targeted clusters (e.g.. iBavo et al.l [201 1[ 
Beia r et al.l 1201 lb . Spectroscopic confirmation of photometric 
candidates in all the regions is needed, however, to construct an 
unbiased IMF, since the contamination by field dwarfs, giants, 
or even galaxies at the faint magnitude ranges being probed is 
known to be high. This is the motivation behind the observations 
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presented in this work. From the initial list of 1 10 brown dwarf 
candida te members of p Oph uncove red through the WIRCam 
survey dAlves de Oliveira et al J 120101) . we had confirmed spec- 
troscopically seven new members in our previous pilot study, six 
of which were classified as brown dwarfs and one very low mass 
star. Here, we present the results of an extended spectroscopic 
follow-up of 34 additional targets. With this survey, we have ob- 
served spectroscopically 70% of our candidates that fall within 
the photometric completeness limit of the WIRCam survey, and 
are nearly complete (~93% of candidates observed) down to a 
few Jupiter masses, through a visual extinction of 20 mag. 

In Sect. |2] we describe the new spectroscopic observations 
and data reduction procedure. Section [3] explains the method 
adopted for deriving the spectral type and extinction for each 
candidate member, and the final results are presented in Sect. [4] 
We discuss the properties of the newly found members and their 
comparison to the cluster's population through Sects. [5] and [6] 
and present the final conclusions in Sect. [7] 

2. Observations and data reduction 

2.1. The ISAAC I VLT spectroscopic follow-up 

We have conducted a near-infrared (near-IR) spectroscopic 
follow-up of 34 targets from our WIRCam/CFHT list of can- 
didate members of p Oph, randomly selected with magnitudes 
brighter than the WIRCam survey completeness limits at 7=20.5 
and //=18.9 mag, and through a visual extinction of approxi- 
mately 20 magnitudes (see Sect. I6.2.TI and Pig . |9j . We addition- 
ally observed five previously known brown dwarfs with spectral 
types determined in the optical, which we use to test our anal- 
ysis methods (see Tables Q] and The observations were car- 
ried out in visitor mode from 4 to 6 May 2 010 with ISAAC 
(Infra red Spectrometer And Array Camera, iMoorwood et alJ 
Il998i) mounted on the Unit Telescope 3 (UT3) of the VLT, 
Paranal (Table Q]). Each target was observed with the low- 
resolution grating in the short wavelength arm of ISAAC and 
a long slit with a width of 0.6" or 0.8" to better match the see- 
ing conditions. Given the faint J magnitudes of most of the tar- 
gets, observations were done only in the H (1.42-1.82/im) and 
K ( 1. 82-2. 5//m) bands. Individual exposure times were scaled to 
the target's brightness and night conditions, and repeated in an 
ABBA pattern to allow for sky subtraction. The slit position was 
aligned with the parallactic angle, except for targets CFHTWIR- 
Oph 20 and 98, which were aligned with a nearby faint target. 
Standard AO stars were observed at regular intervals during the 
night, disk 

The data was reduced with IRAF (Image Reduction and 
Analysis Facility) for all targets by following standard proce- 
dures. The 2D spectra were flat-fielded and corrected for bad 
pixels, after which the sky subtraction was done by subtracting 
adjacent pairs of spectra. The spectra were corrected for dis- 
tortions along the spectral and spatial directions, wavelength- 
calibrated, and then aligned and co-added. Individual frames 
where the spectrum was strongly affected by bad pixels or cos- 
mic rays were removed before the co-adding operation. The ex- 
posure times per individual integration and the final number of 
combined frames for each target are given in TableQ] Each spec- 
trum was extracted with the APALL routine in IRAF. Each ex- 
tracted spectrum was further divided by that of the standard AO 
stars observed at similar airmass to remove the telluric absorp- 
tions, and the relative flux recove red by multipl ying it by a the- 
oretical spectrum of an AO star dPicklesl [l998, taken from the 
ESO webpage). A linear interpolation was made to remove the 



strong intrinsic hydrogen absorption lines from the spectra of the 
AO standard stars, across the lines that are more predominant at 
this resolution (the Brackett series lines at 1.54, 1.56, 1.57, 1.59, 
1.61, 1.64, 1.68, 1.74, 2.17 pm). To calibrate the relative fluxes 
of the spectra originating from the two grisms, we calculated 
the pseudo spectral flux for each band by convolving the spectra 
with the respective WIRCam transmission curves and integrat- 
ing the flux, which are then compared to the flux ratios for H and 
K s band measurements of the WIRCam photometry, and scaled 
accordingly. 

2.2. The Osiris I GranTeCan observations 

We obtained optical spectra for three targets from the WIRCam 
survey, two of them were also observed in the near-IR during the 
ISAAC/VLT run previously described (CFHTWIR-Oph 30, 31), 
while the other had b een observed with Sofl/NTT a nd presented 
in our survey paper yesde Oliveir a et al. IHM CFHTWIR- 
Oph 96). The targets were observed with OSIRIS (Cepa et al. 
2000), a spectrograph operating from 0.365 to 1 .05 /mi, mounted 
on the Gran Telescopio Canarias (GTC), the 10.4 m telescope 
at the Observatory del Roque de Los Muchachos (La Palma, 
Canary Islands). We obtained low-resolution spectroscopy for 
all targets with the R300R grism and an 1" slit aligned with 
the parallactic angle. The observations were obtained in ser- 
vice mode on 22 and 23 May 2010, as well as the spectro- 
photometric standards, flats, bias, and arcs necessary for the cal- 
ibrations. Each target was observed with two exposures of 450 s 
each for CFHTWIR-Oph 31 and 96, and 1500 s for CFHTWIR- 
Oph 30. The data was reduced using standard IRAF routines, 
including bias subtraction, flat-field correction, extraction of the 
ID spectrum, and wavelength calibration. The instrumental re- 
sponse was corrected using spectrophotometric standards, and 
the individual spectra were median-combined. 

3. Methods 

3.1. Spectral classification in the near-IR 

It has been shown by several studies that low-resolution near- 
IR spectroscopy can be used to accurately determine spectral 
types of young stellar objects and confi rm their low gravity (e.g., 
Alter s et alJl2007t lWeights et al. 2009). One of the most evident 
youth indicators is the water vapour abs orption bands that result 
in the triangular shape of the //-band dLuhman & Riekd 1 1 999; 
iLucas etalJl200lt iMcGovern et al.ll2004l) . with a peak redwards 
of that of field dwarfs that instead show a plateau at the same 
wavelength range. The same effect of H2O absorption is seen 
in the K band (>2.3//m), where CO absorption is also present. 
Provided the spectrum has a reasonable signal-to-noise ratio 
(S/N~30), the Na I feature at 2.2//m can be used as a youth in- 
dicator, since it is less prominent in young objects than in field 
dwarfs. Finally, near-IR emission lines characteristic of outflows 
(e.g., H2 1-0 S(0) emission at 2.12 pm) or accretion (e.g., hy- 
drogen recombination line Bry at 2.16 /mi) can be detected even 
at low resolution and provide further evidence of the member- 
ship and youth of a candidate. In this section, we describe the 
methodology developed to investigate the nature of the candi- 
date members of the p Ophiuchi cluster. 

3.1 .1 . Empirical grid of template spectra 

To derive the spectral type and extinction of the candidate mem- 
bers observed with ISAAC, we compare each spectrum to an em- 
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Table 1. Journal of the observations. 



Date 


Slit (") 


CFHTWIR-Oph target 


Known brown dwarfs from literature 


ISAAC/VUF 


2010 May 04 


0.8 


7 (120^x5, 1205X10) 
9 (605X9,605X10) 
12(605X9,605X10) 

15 (405X2, 205X2) 

16 (405X2, 305X2) 

17 (605X8, 605X4) 
31 (205X2, 105X2) 
33 (1205X9, 1205X9) 
66 (605X10, 605X8) 
77 (605X9, 605X10) 
90 (605X8, 605X4) 
100 (605X 8, 605X9) 
103 (605X9, 605X10) 


[L2007b] ChaJl 1070768-7626326 (120.5X14, 1205X10) 
OTS44 (605X10, 605X 8) 


2010 May 05 


0.6 


52 (605X2, 605X2) 
98 (605X15,605X12) 






0.8 


19(120sxl2, 1205X10) 
20(1205X9, 1205X10) 
30 (405X2, 205X2) 
37 (405X4, 205X4) 
42(1205X16, 1205X14) 
56 (605X2, 605X2) 
69 (205X2, 205X2) 
78 (605X4, 305X4) 
107 (205X2, 205X2) 


CHSM17173 (605X6, 605x6) 
[NC98] Cha HA 11 (605x7, 605x8) 


2010 May 06 


0.6 


1 (1205X14, 1205X7) 
3 (605X2, 605X2) 
6 (1205X8, 605X5) 
18 (1205X13, 1205X14) 
22 (605X4, 605X4) 
28 (1205X20, 605X10) 
53 (605X7, 605X4) 
101 (1205X6, 605X6) 
104 (1205X6, 605X6) 






0.8 


58 (605X2, 605X4) 


TWA 26 (605X 8, 6O5X8) 


Osiris/GTC 


2010 May 22 


1.0 


31 (4505X2) 




2010 May 23 


1.0 


30(15005X2) 
96 (4505X2) 





Exposure times are given in seconds per individual integration, times the number of total integrations for the H and K grisms, respectively. 



pirical grid of near-IR template spectra of young stellar objects 
while varying the extinction. The objects that integrate the grid 
were chosen out of a list of the spectra of 37 confirmed members 
of Taurus and IC 348 (l-3Myr), which have an approximate age 
to that of p Oph, and have spectral types from M3 to M9.5 clas- 
sified in th e optical through comparison to averages of dwarfs 
and gia nts dBriceno et al.ll2002tlLuhman et al . 2Q03a'S lLuhrrianl 
l2004bl) . To determine the extinction towards each of the tem- 
plate spectra, the spectro-photometric colours were calculated 
by converting to magnitudes the integrated flux of each spectra 
convolved with the 2MASS transmissio n filter curves, and using 
the respective fluxes for zero magnitude dCohen et al.l2 003). The 
extinction (Ay) was determined by comparing t hose to the ex- 
pecte d J-H and H-K colours for young objects (Luhma n et alJ 
2010, see Table 13), and applied to deredden each spectrum ac- 
cordingly. 

Based on their photometric properties and evolutionary mod- 
els, the candidate members we observed could have spectral 
types as late as L5. Therefore, we added to this list near-IR 
spectra (K. Luhman, priv ate communicatio n) of optical stan- 
dard young-field L dwarfs dCruz et aTl f2009) with spectral types 
from L0 to L5 and spectral features indicative of low gravity. 



Although the membership of these objects to nearby loose as- 
sociations has yet to be confirmed, the low-gravity features in 
their spectra are consistent with an age of *T0Myr. Even tak- 
ing the most conservative approach to this study, the behaviour 
of the gravity-sensitive indicators in their spectra (mainly from 
the alkali lines and the VO bands, see Kirkpatrick 2008) im- 
ply an upper-limit for their age of lOOMyr making them more 
suitable for spectral classification of young brown dwarfs than 
old field dwarfs. Table [2] lists the objects selected to build the 
grid of young templates, including their identification, together 
with the extinction values derived for the members of IC 348 and 
Taurus. 2MASS photometric errors (between 0.05 to 0.11 mag) 
are approximately the same as the differences between photo- 
metric and synthetic colours. Intermediate subclasses have been 
calculated by averaging between subsequent spectral types. The 
second dimension of the fitting matrix is created by progres- 
sively reddening all templates in the grid b y steps of 0.1 visual 
magnitude, using the extinction law of lFitzpatrickl (fT999). 



An analogous grid was constructed that is composed of 
field (1-5 Gyr) template spectra of G and K stars dRavner et alJ 
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2009) retrieved from the IRTF Spect ral Libraryfl, together 
with M and L dwarf optical standard s (Bureasser et al. 2004; 
Burgasser & McElwain 2006; Burgasser et al. 2008) retrieved 
from the SpeX Pris m Spectral Librarie all taken with the 
SpeX spectrograph dRavner et alj 12003b mounted on the 3 m 
NASA Infrared Telescope Facility. These are reddened in the 
same way as described for the grid of young objects, and the grid 
is smoothed by calculating intermediate spectral types from the 
averages of the templates. This secondary grid is used to double- 
check a posteriori that the youth indicators described above are 
present in the candidate spectra. 



3.1.2. Fitting technique 

For each candidate spectrum, we determine the best fit in spec- 
tral type and extinction by minimizing a go odness-of-fit statisti c 
between the candidate and each template dCushing et al.l l2008). 
defined as 



(f: 



■candidate 



-^template \ 



andidate 



(i) 



where is the w eight of each spectr al element (taken as its 
spectral width as in iCushing et al.ll2008l) . fx and are the flux 
densities of the candidate and template, respectively, and a\ 
is the noise for the candidate spectrum. Regions of the spec- 
tra dominated by telluric absorption were excluded from the 
fit, which was done over the wavelength ranges 1.45-1.8 and 
2-2.45 ftm. Since the ISAAC spectra cover only the H and K 
bands, we have used the photometric point of the J band in the 
fit, which is compared to the spectro-photometric magnitude de- 
rived by each template in the grid, and given a weight that is 
the sum of the individual weight values in the J band of the 
template. This approach is necessary in order to minimise the 
effect the degeneracy between extinction and spectral type have 
on the overall slop e of a young late-M and L-type spectrum (e.g., 
lAllers et al.ll2007l) . 

3.1.3. Testing the fitting method 

We tested this technique by determining the near-IR spectral 
type and extinction for the 30 aforementioned confirmed young 
members of Taurus and IC 348, where the /-band part of each 
spectrum is replaced by its integrated flux to mimic our ISAAC 
dataset. The vast majority of the derived spectral types agree to 
better than 1 subclass with those determined from the optical, 
and to better than 1 magnitude of Ay from the extinction values 
determined from the near-IR colour-colour diagram, as shown 
in the histogram of the parameter's deviation in Fig.Q] We take 
these to be the intrinsic errors in the fitting method and to con- 
firm that an accurate spectral type can be derived when replac- 
ing the information of the J band spectra with its integrated flux. 
Therefore, the final uncertainties in the spectral type and extinc- 
tion derived for the p Oph candidate members observed with 
ISAAC/VLT are taken to be 1 spectral subclass and 1 mag. 

We repeated the test for the known brown dwarfs in the liter- 
ature observed during our ISAAC run. Four out of the five brown 
dwarfs we observed have a spectral type solution from our fitting 
technique that agrees to within one subclass in respect to their 
published optical spectral type. For the brown dwarf Cha Ho- 1 1 
we find a spectral type that is ~2 subclasses earlier. The derived 



CD 




-4-3-2-101234 
Deviation 

Fig. 1. Histogram of the deviation between the near-IR spectral 
type (solid line) and the extinction (dotted line) derived with our 
fitting technique for young members of Taurus and IC 348, and 
the spectral type derived from the optical and the extinction de- 
rived from the near-IR colour-colour magnitude diagram. 

values of extinction show a dispersion of <2.6 mag to those pub- 
lished, which are zero for all targets. This could be explained by 
our reliance on 2MASS photometry for the relative calibration 
of the spectral bands, propagating the photometric errors to the 
spectrum. Furthermore, particular geometries of disks around 
the brown dwarfs could cause variability at near-IR wavelengths, 
as well as contribute to the K-bwd overall flux, which would im- 
pact the extinction derived, but not significantly affect the overall 
spectral features. Table[3]summarises the properties of these ob- 
jects as retrieved from the literature, as well as the parameters 
derived from our spectral fitting method. 

3.2. Spectral classification in the optical 

The most distinct feature in the optical spectra of M dwarfs are 
the absorption bands from titanium oxide (TiO) and vanadium 
oxide (VO), which are very sensitive to temperature, increasing 
in strength as we move to later spectra l types in a smooth pro- 
gression (e.g., Kirkpatri ck et alj |l993). To determine the spec- 
tral type and extinction of our targets, we employ the same algo- 
rithm as is described in the previous section, using the goodness- 
of-fit statistic to determine the best match between the candi- 
date member and the grid of templates. The grid is composed 
of young objec ts that are members of Chamaeleon (~lMyr, 
Luhman 2004a), which have been shown to be suitable for rel- 
ative spectral classificat ion of young M dw arfs in star-forming 
regions (see Table 2 in Rid dick et al.l l2007). The templates are 
reddened progressivel y in steps of 0.1 v isual magnitude, using 
the extinction law of iFitzpatrickl d!999l) . The typical errors of 
this method in the spectral type domain we are studying (late M 
to early L) are of the order of 1 subclass in spectral type and 
1 mag nitude of extinction (e.g jGuieu et al.f2006l:lRiddick et alj 



2007). 



4. Results 

The numerical spectral fitting technique has been applied to all 
yt- andidate members observed with ISAAC/VLT and Osiris/GTC, 
and the solution for spectral type and extinction is summarised 



1 Available at http://irtfweb.ifa.hawaii.edu/~spex/IRTF_Spectral_Library)t andidate members observed with ISAAC/VLT and Osiris/GTC, 

2 Available at http://www.browndwarfs.org/spexprism/ 
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Table 2. Empirical grid of young template spectra. 



2MASS Designation 


Region 


Optical SpT 


Ref." 


Ay" (mag) 


(J-H)2MASS C 


{H-K S )2MASS C 


(J-H)wiRCam d 


(H-K s )wiRCam d 


J03442812+3216002 


IC348 


M3.25 


1 


2.0 


0.69 


0.22 


0.70 


0.20 


J04311578+1820072 


Taurus 


M4.25 


2 


0.4 


0.63 


0.26 


0.63 


0.24 


J04 183030+2743208 


Taurus 


M5.5 


2 


0.6 


0.56 


0.29 


0.54 


0.27 


J04 180796+2826036 


Taurus 


M6 


3,4,5 


1.6 


0.55 


0.32 


0.50 


0.31 


J04394748+2601407 


Taurus 


M7 


6,2 


5.4 


0.60 


0.38 


0.55 


0.37 


J043057 18+2556394 


Taurus 


M8.25 


2 


0.7 


0.66 


0.46 


0.60 


0.45 


J04190126+2802487 


Taurus 


M9 


7 


1.3 


0.74 


0.56 


0.65 


0.55 


J01415823-4633574 


field 


L0 


8 


0.0 


1.01 


0.65 


0.91 


0.65 


J22081363+2921215 


field 


L3 


8 


0.0 


1.19 


0.78 


1.04 


0.78 


J05012406-0010452 


field 


L4 


8 


0.0 


1.31 


0.82 


1.12 


0.82 


J03552337+ 1133437 


field 


L5 


8 


0.0 


1.57 


1.07 


1.42 


1.08 



(fl> S pe ctral types from: LlLuhman etall d2003bl). 2:lBriceno et all d2002h. 3 lStrom & Stroml d!994h . 4 lLuhman et"a!] ( ll998h , 5 lBriceno et all 
d!998h . 6 lMartm et all d200lh . 7: lLuhman et alj j2003J l. 8: ICruz et ail d2009h . (b) Visual extinction determined from the J-H vs. H-K s diagram, 
except for the field L dwarfs. w Colours computed from the dereddened spectra in the 2MASS photometric system. <£/) Colours computed from 

the dereddened spectra in the WIRCam/CFHT Vega system. 



Table 3. Brown dwarfs in the literature. 



2MASS Designation 


Region 


Opt SpT" 


Ref. fl 


Other name 


NIR SpT" 


Ay" (mag) 


Jl 1082927-7739198 


Chamaeleon I 


M7.25 


1 


[NC98] Cha HA 1 1 


M5.5 


1.5 


Jl 1102226-7625 138 


Chamaeleon I 


M8 


1 


CHSM 17173 


M7.5 


0.8 


J11395113-3159214 


Twa Hydra 


M8 


2 


TWA 26 


M7.25 


1.8 


Jl 1100934-7632178 


Chamaeleon I 


M9.5 


3 


OTS44 


M9.25 


2.6 




Chamaeleon I 


L0 


4 


[L2007b] Cha Jl 1070768-7626326 


M9 


0.8 



(a) Spectral types from: 1: Luhman (2004a), 2: lRice et" al] d2010l) , 3: lLuhmanet~aHd2004h , 4: lLuhman etaU d2008al) . (w Near-IR spectral type 

and extinction derived from numerical fitting. 



in Table|4] The final dereddened spectra are shown in Figs. [21 El 
and|4] In the following sections we discuss the fits to the candi- 
dates in more detail, including the membership assessment and 
the level of contamination in the observed sample. 

4.1. New spectroscopically confirmed members of p Oph 

By applying the numerical spectral fitting technique to all spec- 
tra, we classify 20 of the observed targets as new young mem- 
bers of the p Ophiuchi cluster, with spectral types between M3 
and L4. An additional target, CFHTWIR-Oph 96 had been pre- 
viously classified as a member through near-IR spe ctroscopy 
(M8.25, A v =1.10 mag. lAlves de Oliveira etal] I2010I) . In this 
work, we derive a more accurate optical spectral type (M7.75, 
A(/=2.4 mag), which agrees, within the errors, with the near-IR 
spectral type. We classify 14 targets as likely contaminants. 

4.1.1. Membership 

The youth of each target, hence membership, was assigned based 
on the presence of signatures of low-gravity or/and accretion 
disks, as well as each target's spatial location in relation to the 
cloud environment. For that end, the results from the spectro- 
scopic follow-up were complemented with the available infor- 
mation in the literature, most of which had already been de- 
scribed and su mmarised in Table 1 of the W IRCam survey pa- 
per of p Oph dAlves de Oliveira eta fl l2010l) . For targets with a 
near-IR spectral type later than ~M3, the shape of the water- 
vapour absorption bands has been matched to that of spectra of 
young stellar objects with low gravity, providing a clear indica- 
tion of their youth (denoted as H2O in Table [4]). The detection 
of mid-IR excess in the Spitzer photometry consistent with that 



expected for young objects surrounded by disks (ex in Table |4| 
is another indication of membership. For these objects, accre- 
tion might still be on-going, so any signatures of Br y emission 
are also indicative of youth, as well as H2 emission arising from 
outflows. Furthermore, all objects selected as members do not 
show any sign of Na I absorption, which is consistent with them 
having low gravity. Finally, the spatial location of an object also 
contains relevant information of whether it belongs to the cloud 
or not. For targets that do not show strong IR excess, the extinc- 
tion derived in the numerical spectral fitting is not contaminated 
by the contribution of a disk, and should reflect the position of 
the target in respect to the depth of the cloud. We compared that 
value to the extinction map of the Ophiuchus cloud provided 
by the COMPLETER project dRidee et al.l|2006h . computed with 
the NICER algorithm dLombardi & Alvesl2001l) using data from 
the Spitzer Lega cy Program "From Molecular Cores to Planet 
Forming Disks" (Evans et al. 2003). The extinction provided by 
the map is determined over a beam that is larger than our tar- 
get observations and also over the entire depth of the cloud. 
Nonetheless, it still provides a good estimate of the environ- 
ment for each candidate. All of our candidate members (except 
CFHTWIR-Oph 101, see discussion below) have extinctions that 
place them inside the cloud, when compared to the Spitzer ex- 
tinction map. That is not the case, for example, for the M dwarfs 
we listed in Table [4] as likely contaminants, which have extinc- 
tions equal to or even greater than the expected values for the 
respective location placing them behind the cloud. 

From the targets classified as members, sixteen (CFHTWIR- 
Oph 16, 18, 30, 31, 37, 53, 56, 58, 66, 78, 90, 96, 98, 103, 
104, 107) show mid-IR excess indicative of a disk, as well as 
spectral features consistent with low gravity, which we con- 



3 Available at http://www.cfa.harvard.edu/COMPLETE 
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Table 4. Spectral type and Ay determined through numerical spectral fitting. 



CFHTWIR-Oph" 


RA 


Dec 


J (mag) 


H (mag) 


K s (mag) 


SpT" 


Ay (mag) 6 


T(K) C 


Properties 


Confirmed new members 


9 


16:26:03.28 


-24:30:25.9 


17.76 


16.33 


15.32 


L0 


6.6 


2250 


H 2 0, CO 


16 


16:26:18.58 


-24:29:51.6 


17.09 


14.92 


13.55 


M6.75 


18.8 


2908 


ex, H 2 0, vciv 


18 


16:26:19.41 


-24:27:43.9 


18.92 


17.16 


16.04 


L0 


7.8 


2250 


ex, H 2 0, CO 


30 


16:26:36.82 


-24:19:00.3 


16.24 


14.35 


13.06 


M3.25 


4.9 


3379 


ex, var, X-ray 


31 


16:26:37.81 


-24:39:03.3 


14.65 


13.45 


12.66 


M5.25,M5.5 


5.7, 8.9 


3091 


ex, VO 


33 


16:26:39.69 


-24:22:06.2 


18.16 


16.74 


15.68 


L4 


3.6 


1650 


H 2 0, CO 


37 


16:26:40.84 


-24:30:51.1 


17.32 


14.77 


13.18 


M5.75 (M7.5) rf 


23.4 (22.) 


3024 


ex,H 2 


66 


16:27:14.34 


-24:31:32.0 


18.42 


16.53 


15.30 


M7.75 (M4.5) d 


15.1 (17.6) 


2753 


ex, H 2 0, var 


77 


16:27:25.64 


-24:37:28.6 


18.23 


16.54 


15.36 


M9.75 (M7) rf 


10.0(15.4) 


2263 


H 2 0, CO 


78 


16:27:26.23 


-24:19:23.1 


16.35 


14.36 


13.04 


M7.75 


16.4 


2753 


ex, H 2 0, var 


90 


16:27:36.59 


-24:51:36.1 


16.83 


15.65 


14.85 


L0 


2.4 


2250 


ex, H 2 0, CO 


96 


16:27:40.84 


-24:29:00.8 


14.60 


13.76 


13.19 


M7.75 


2.4 


2753 


ex, VO 


98 


16:27:44.20 


-23:58:52.4 


17.08 


15.78 


14.98 


M9.75 


3.0 


2263 


ex, H 2 0, CO 


100 


16:27:46.54 


-24:05:59.2 


17.94 


16.34 


15.26 


L0 


8.0 


2250 


H 2 0, CO 


101 


16:27:47.25 


-24:46:45.9 


19.26 


16.48 


14.78 


M7 


25.8 


2880 


H 2 0, em 


103 


16:28:10.46 


-24:24:20.4 


17.74 


16.15 


15.07 


L0 


7.5 


2250 


ex, H 2 0, CO 


107 


16:28:48.71 


-24:26:31.8 


14.31 


13.60 


13.13 


M6.25 


2.3 


2963 


ex 


Members with uncertain spectral type 


53 


16:26:56.36 


-24:41:20.6 


18.60 


16.34 


14.92 


L3? 


8.6 




ex, H 2 0, CO 


56 


16:27:03.59 


-24:20:05.6 


17.13 


15.01 


13.72 


M3? 


20.1 




ex, var, em 


58 


16:27:05.94 


-24:18:40.3 


17.05 


15.77 


14.85 


L3? 


1.5 




ex, H 2 


104 


16:28:11.61 


-24:37:29.9 


19.51 


16.95 


15.41 


L3? 


10.5 




ex, H 2 


Likely contaminants 


1 


16:25:19.15 


-24:19:27.4 


19.07 


17.33 


15.86 








ex 


3 


16:25:23.78 


-24:22:58.4 


16.30 


14.85 


13.97 


M5 


9.7 




Nal 


6 


16:25:46.63 


-24:23:36.4 


18.87 


16.69 


15.12 








ex 


7 


16:25:46.98 


-24:11:54.1 


18.34 


16.85 


15.68 








ex 


12 


16:26:09.99 


-24:14:40.1 


18.15 


16.35 


15.25 


M4 


15.3 




Nal 


15 


16:26:16.33 


-24:39:30.7 


15.79 


14.15 


13.11 


M6.5 


12.4 




Nal 


17 


16:26:19.21 


-24:41:30.1 


17.44 


15.80 


14.61 










19 


16:26:22.27 


-24:37:09.5 


19.15 


17.10 


15.86 










20 


16:26:22.89 


-23:58:34.2 


18.01 


16.73 


15.80 










22 


16:26:25.10 


-24:41:32.7 


17.03 


15.52 


14.54 










28 


16:26:35.28 


-24:42:39.2 


19.68 


17.79 


16.36 








ex 


42 


16:26:44.51 


-24:04:08.8 


18.64 


17.19 


16.08 










52 


16:26:55.54 


-23:57:36.5 


16.02 


15.08 


14.39 








ex 


69 


16:27:18.53 


-24:07:22.0 


14.76 


13.55 


12.79 


M7.5 


6.7 




Nal 



(fl) Name as in lAlves d e Oliveira et al.l d2010t) . where additional information regarding each target is included. {b) Based on the tests performed to 
the spectral fitting method, uncertainties on the spectral type and A v are estimated as 1 spectral subclass and 1 mag, respectively. (c) See Sect. 16.21 
where the temperature scale applied is explained. {d) See Sect. l4.1.2l for the discussion. 



sider the most constraining clues to their membership. In ad- 
dition, five from those (CFHTWIR-Oph 16, 30, 56, 66, 78) 
have the near-IR variability typ ical of young stellar objects 
dAl ves de Oliveira & Ca sali 2008). The membership of the five 
targets without mid-IR excess (CFHTWIR-Oph 9, 33, 77, 100, 
101) is assigned based on complementary diagnostics and dis- 
cussed below in more detail. 

CFHTWIR-Oph 9, 33, 77: Candidate 9 is located in a region of 
the cloud with high extinction (Ay* 17 mag), and the derived ex- 
tinction value for this object is Ay^6.6 mag, clearly placing it in 
the middle of the cloud. Its K band spectra also clearly shows the 
CO overtones characteristic in spectra of early L dwarfs, which 
is consistent with its near-IR and optical photometric colours. 
CFHTWIR-Oph 33 is located close to the core of the cloud in 
a region with an Ay*24 mag, this target shows very peaked H 
and K band spectra characteristic of young objects, and with an 
Ay^3.6 mag it is definitely not a background source. It is too 
bright to be a late L foreground field dwarf that its spectra in- 
dicate, further confirming its membership to the cluster. Similar 



to these two targets, CFHTWIR-Oph 77 is located in an area 
of high extinction (Ay~25 mag), and following the same argu- 
mentation, both the derived Ay and spectral type are consistent 
with it being a member of the cloud. It could not be reproduced 
by any of the templates in the grid of field dwarfs. The spectral 
classification of this target is discussed further in Sect. 14.1.21 

CFHTWIR-Oph 1 00 and 101: The target 100 is in a region of 
relatively low extinction (Ay^5-8 mag), which places the object 
at the further edge of the cloud, when compared to the extinction 
derived by numerical fitting (Ai/»8 mag). When derredening the 
colours of the target by this amount, the object still lies confi- 
dently above the main-sequence and is rather consistent with the 
expected colours of an early L young brown dwarf. The signal- 
to-noise ratio in the K band is not enough to make a definite as- 
sessment of the Na I feature. The CO overtone absorptions in the 
K band further support its late spectral type, which again would 
be too bright to be a field dwarf behind the cloud. CFHTWIR- 
Oph 101 is the only target we classified as a member for which 
the derived extinction value exceeds that expected for the region 
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CFHTWIR-Oph 31 (M5.25) 
I 

CFHTWIR-Oph 37 (M5.75) 



CFHTWIR-Oph 107 (M6.25) 
CFHTWIR-Oph 16(M6.75) 



CFHTWIR-Oph 101 (M7) 
Oph 66 (M7.7 



CFHTWIR-Oph 66 (M7.75) 



CFHTWIR-Oph 78 (M7.75) 



'*~r> 



I I I I 1 I I I 1 I I I 1 I I 1 

CFHTWIR-Oph 98 (M9.75) 


1 1 1 1 1 1 1 1 1 1 


CFHTWIR-Oph 77 (M9.75) 


- 


■ " J\ 

V 

CFHTWIR-Oph 9 (L0) 




/A 

" CFHTWIR-Oph 18 (L0) 




CFHTWIR-Oph 90 (L0) 




* /\ 

CFHTWIR-Oph 100 (L0) 




" " J\ 

CFHTWIR-Oph 103 (L0) 




FT 1 
" CFHTWIR-Oph 33 (L4) 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 



1.2 1.4 1.6 1.8 2.0 2.2 2.4 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

A,(jim) 

Fig. 2. IS AAC/VLT low-resolution spectra of the new confirmed 
members in p Oph. All spectra are corrected for extinction with 
the values found through the numerical fitting procedure. 



3.0 

I 2-5 

"cO 

<§ 2.0 
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CO 

I 1.0 
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z 0.5 



CFHTWIR-Oph 56 (early M) 




CFHTWIR-Oph 53 (late M/early l_T 



1.2 1.4 1.6 1.8 2.0 2.2 2.4 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

A.(jim) 

Fig. 3. IS AAC/VLT low -resolution spectra of the members with 
uncertain spectral type. 



from the extinction map. When looking at the near-IR images, 
it is, however, clear that differential extinction is present, since 
this part of the cloud shows a ridge-like structure. The target sits 
on one of the dense ridges, which seems to explain the high ex- 
tinction, and has most likely been smoothed in the low spatial 
resolution of the extinction map. There is no indication of Na I 
absorption, further supporting its low-gravity, while there is an 
indication of Br y emission (2.16 pm) showing that accretion 
may be on-going. Finally, the peaked shape of the H band is 
not reproduced by any spectra of the grid of field dwarfs. Based 



1 1 1 1 I 1 1 1 1 I 

CFHTWIR-Oph 30 (M3.25) 



I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 




Fig. 4. Osiris/GTC optical spectra of three confirmed members 
in p Oph. All spectra are corrected for extinction with the values 
found through the numerical fitting procedure. The best-fitting 
template is overplotted in red. 



on these arguments, we classify both targets as members of the 
p Ophiuchi cluster. 

4.1.2. Comments on individual sources 

CFHTWIR-Oph 30: The candidate member CFHTWIR- 
Oph 30 has been previously associated with the cloud by several 
diagnostics, including mid-IR excess emission consistent with 
that of YSOs (e.g., iBontemps et al.l 12001b iGutermuth et al.l 
1 20091) and X-ray emission (iGagne et al.1 12004) . More recently, 
lAlves de Oliveira & Casalil (|2008) found this object to show 
extreme variability in the near-IR, with a brightening of 3.5 
magnitudes in the K band in only one year, which the authors 
suggest could indicate an EXor-type object, a class of young 
variable stars where the large magnitude variations are thought 
to be caused by massi ve infall of circumstellar material onto 
the central star (e.g., [Herbig 2008). Based on the acquisition 
images of the ISAAC dataset, we have estimated a magnitude 
of ZT=14.54±0.3, which would indicate that in 2010 the target 
was again approaching its minimum state. As a result, the 
necessary spectroscopic integration times have been underes- 
timated and the signal-to-noise of the ISAAC spectra is too 
low for a reliable spectral fitting. This is underpinned further 
by the large uncertainty of the near-IR magnitudes at the time 
of the spectroscopic observations, crucial for calibrating the 
relative spectral fluxes of the various bands. This target was also 
observed spectroscopically in the near- IR bv lGeers et all d201 lb 
(source #3 in their Table 1) and classified as non-substellar. We 
obtained an optical spectrum of CFHTWIR-Oph 30 with GTC, 
allowing more reliable analysis. We derive a spectral type of 
M3.5 with a visual extinction of 4.7 magnitudes (see Figj4]i. 
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There is no indication of Her emission in the optical spectrum, 
which could support the scenario where the YSO is at its rest 
accretion state. 



We have compiled the available IR 
from the literatur e, which are described 
articl es: iBarsonv et al.l 
(I2003L 



lowing 



ICutri et alJ (I2QQ3L 2MASS} 1 
d2008[ WFCAM/ UKIRT) 



lfl997i 



measurements 
in the fol- 
SQIID/KPNO) 



Alves de Oliveira & C asali 
Alves de Oliveira et all (1201 Oi 



WIR Cam/CFHT ) iGeersetalJ (1201 lL MOIRCS/Su baru) 



iGutermufh etail (I2009L IRAC/Spitzer). and I Wright etaD (I20i6L 
WISE). The resulting light curve for the different bands is 
shown in Fig. [5] Even though we did not correct the mag- 
nitudes to account for the differences between photometric 
systems, this cannot account for the large amplitude that is 
seen, and it can only be explained by intrinsic variability. This 
plot strengthens the rare c olour effect already measured by 
lAlves de Olive ira & Casalil (2008), with the star bluer when 
fainter and the cause of variation possibly related to the star's 
disk. For example, as a clump of material or warp in the disk 
moves out of the observers line-of-sight, the star-disk system 
becomes fainter in the IR since less light is being absorbed and 
reprocessed by the circumstellar disk. Disks radiate strongly 
at longer wavelengths, so when the warped inner disk region 
rotates out of the observed fraction of the star, the stellar colours 
become bluer, approaching its intrinsic colours. However, the 
intrinsic fluxes of an M3.5 YSO are much brighter than what is 
observed at its minimum. Indeed, its optical spectral type would 
indicate that the nominal state of the target is at its brightest 
stage. If taking this into account, the variability could instead 
be explained by a companion object, since the percentage of 
coverage of a cold spot on the surface needed to reproduce 
such behaviour is not physical when considering the magnitude 
variations. Then, the observed behaviour could be due to a 
system consisting of two objects with similar spectral type (i.e. 
similar intrinsic luminosity), one having strong IR excess and 
undergoing variable circumstellar extinction by the disk (i.e., 
close to edge-on), the other component either disk free or at 
least not obscured by its disk (e.g. different inclination). Then, 
the whole system would become fainter in the near-IR when the 
disk obscures the primary. At that time the secondary would 
start to dominate the near-IR integrated light. If the secondary 
has no (or a smaller) IR excess, its near-IR colours will be bluer 
than that of the primary, and the system would in turn become 
bluer when fainter in the near-IR. A detailed modelling of the 
system is beyond the scope of this paper, and will most likely be 
addressed by dedicated monit oring surveys such as YSOVAIiQ 
dMorales-Calderon et alj|201 ll) . 



CFHTWIR-Oph 37, 66, 77: For all the targets we analysed, the 
near-IR numerical fit is done both with the J band photomet- 
ric point and without it, for which we obtain the same solution, 
confirming that the spectral shape of the H and K bands is the 
dominant component of the fit, as opposed to the JHK colours 
on which the spectra are calibrated. This is, however, not the 
case for three targets (CFHTWIR-Oph 37, 66, 77) where we ob- 
tained two different solutions depending on whether or not we 
included the comparison between the J photometric point of the 
target to the integrated flux of the template's J spectra in the 
fit. The different solutions are in fact correlated, since they are 
extreme cases of the degeneracy between extincti on and spec- 
tral type, which exists at this low resolution (e.g., lAllers et"al] 




2449000 2450000 2451000 2452000 2453000 2454000 2455000 



Julian Date 
1.5 



2449000 2450000 2451000 2452000 2453000 2454000 2455000 



Fig. 5. Light curve for CFHTWIR-Oph 30 combining all 
the IR measurements available (SQIID/KPNO, 2MASS, 
WFCAM/UKIRT, WIRCam/CFHT, MOIRCS/Subaru, 
IRAC/Spitzer, and WISE), where uncertainties are typically 
<0. 1 mag. The magnitudes have not been corrected to a unique 
photometric system, though this cannot account for the large 
variations that can only be explained by intrinsic variability. 



Available at http://ysovar.ipac.caltech.edu 



120071) . and also at a modest signal-to-noise of the targets. In the 
case of CFHTWIR-Oph 37 and 66, their spectral energy distri- 
bution (SED) shows a strong IR excess up to 24 fim, meaning 
the contribution from the disk emission could affect the determi- 
nation of the extinction. If we investigate their near-IR colours, 
CFHTWIR-Oph 37 is consistent with the earlier spectral type 
solution (M5.75, Ay=23.4 mag), while the faint magnitudes of 
CFHTWIR-Oph 66 are in agree with the later spectral type so- 
lution (M7.75, A(/=15.1 mag). We adopt these classifications 
throughout the rest of the paper. CFHTWIR-Oph77 does not 
show mid-IR excess, and its colours are consistent with a later 
spectral type, which is also supported by the appearance of the 
CO overtone absorption in its K band spectra, and therefore we 
adopt the spectral type solution M9.75 with Ay=10.0 mag. 



CFHTWIR-Oph 53, 56, 58, 104: These targets are assigned a 
spectral type that carries a larger uncertainty than the other can- 
didates, due to both the poor signal-to-noise ratio and strong 
mid-IR excess (see Fig. gj. CFHTWIR-Oph 53, 58, and 104 
show strong water vapour absorption in the H band, as well 
as pronounced CO overtone absorption in the K band, and are 
largely consistent with a late M to early L spectral type. The 
signal-to-noise in the H band is too low to assess the quality of 
the fit, and the solutions presented rely mainly on the K band 
features. CFHTWIR-Oph 56 appears to have an earlier spectral 
type than M3, the hottest template in our grid. They all have 
clear signatures of youth and are certainly members of the clus- 
ter. We include the best spectral solution found in Table [4] but 
caution that the error in spectral type is significantly higher than 
for the other targets. We take a conservative approach and do not 
include them in the analysis of the brown dwarf population in 
Sects. |5]and|6] 
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4.2. Likely contaminants 

We classify fourteen candidates in our sample as contaminants, 
which are most likely background giants or field M dwarfs. All 
these targets (CFHTWIR-Oph 1, 3, 6, 7, 12, 15, 17, 19, 20, 22, 
28, 42, 52, and 69) have a non-negligible amount of extinc- 
tion, so they cannot be foreground stars. For most of the tar- 
gets (CFHTWIR-Oph 1, 6, 7, 17, 19, 20, 22, 28, 42, 52), there 
is no sign of water absorption features in their spectra, imply- 
ing a spectral type that is earlier than M3. However, a young 
object (~lMyr) at the distance of p Oph (~130 pc), ought to 
be much brighter than the faint magnitudes these targets show. 
Therefore, they are most likely giants and AGB stars, also known 
to have mid-IR excess, which are seen through the cloud. The 
targets classified as field M dwarfs (CFHTWIR-Oph 3, 12, 15, 
69), show spectral features that are better matched by our grid 
of template field dwarfs, such as a plateau in the shape of the 
H band, and the Na I absorption feature characteristic of older 
objects. Furthermore, they do not show any other sign of youth, 
such as mid-IR, and all have extinction solutions that place them 
behind the cloud when compared to those predicted by the ex- 
tinction maps. 



members classified with a spectral type of M6 are listed in the 
table of brown dwarfs as members at the substellar limit. To 
be consistent with other studies of star-forming regions and to 
facilitate comparison of properties across different populations, 
we label young objects with spectral types later than M6 as 
brown dwarfs. By this criterion, a t the time of the latest review 
of the cluster ( Wil king et alJ l2008). only 1 1 brown dwarf mem- 
bers were known in Ophiuchus. With the WIR Cam survey and 
the su bsequent spectroscopic follow-up dAlves de Oliveira et alj 
2010, and this work), we have confirmed 19 new brown dwarfs 
(22 if we include CFHTWIR-Oph 53, 58, and 104, which have 
uncertain spectral types between late M and early L), a threefold 
increase in the subste llar population o f p Op h. Based on the op- 
tical classification of Comeron et al. (2010), we do not include 
GY 10 (M5.5) in the substellar census and update the spectral 
type of GY 1 1 . Other rece nt contributions to the census of the 
cluster have been done by iMcClure et al.l d2010l) . where spec- 
tral types were derived up to M6 for severa l members, and th e 
confirmation of one brown dwarf member bv lGeers et alJ(l201 ll) . 



4.3. Effectiveness of the photometric selection 

We combined the results for all the candidate members selected 
in our WIRCam/CFHT survey that have been observed spectro- 
scopically. These include seven confir med members presented 
in our previous spectroscopic follow-up dAlves de Oliveira et al.l 
1201 Oh . where we found four contaminants (these were excluded 
from the list of cand idate members of the cluster). Furthermore, 
IMcClure et aOJ2010h confirm the membership of six of our can- 
didate members. In this study, we have confirmed as members 
of the cluster 20 additional targets. The candidates classified as 
contaminants total 18. In absolute terms, this implies a success 
rate of 65%. If we take all the candidate members within our 
photometric completeness limits into account (see Fig. by 
including 17 targets that have not been observed spectroscopi- 
cally, we derive lower and upper limits of the success rate of the 
WIRCam/CFHT survey of p Oph of 49% and 74%, respectively. 
This agre es with previous estimates in surveys of other clus- 
tgrsje.g.. loiiveira et al]l2009r. IScholz etall 120091: ILodieu et alJ 
1201 lbl) . and reenforces both the need for spectroscopic confir- 
mation of candidate members selected from photometry and the 
limitations of previous IMF studies based only on photometric 
selected samples. 



5. The brown dwarf population in p Ophiuchi 

5.1. An updated census 

Despite being one of the closest and youngest star-forming 
regions where brown dwarfs should be more easily detected, 
the high and variable extinction of the p Ophiuchi cluster has 
meant that its substellar population remained largely unknown. 
We compiled a list of all brown dwarf members of p Oph 
confirmed spectroscopi cally to date, which is an updated 
version of Table 4 in lAlves de Oliveira et alJ d2010h . At the 
age of p Oph (~lMyr), the evolutionary models of the Lyon 
group, combined with th e empirical temperatur e to spectral type 
relation determined by Luh man et al.l d2003bl) . imply that the 
mass limit for hydrogen burning corresponds to a spectral type 
of -M6.25 (Luhman et all l2007l) . Given the limited accuracy in 
near-IR low-resolution spectroscopy, typically of one subclass, 
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Table 5. Spectroscopically confirmed brown dwarfs in p Oph. 



Identifier 


RA 


Dec 


SpT 


Ay (mag) 


References" 


J (mag) 


H (mag) 


K s (mag) 


[3.6] (mag) 


[4.5] (mag) 


[5.6] (mag) 


[8.0] (mag) 


[24] (mag) 


Members at the substellar limit 


GY92 15 


16:26:22.98 


-24:28:46.1 


M6 


10.7 


8 


14.93±0.04 


12.80±0.02 


11.53±0.02 


10.41±0.06 


10.00±0.06 


9.50±0.07 


8.46±0.09 


6. 11 ±0.63 


GY92 109 


16:26:42.89 


-24:22:59.1 


M6 


13.5 


8 


15.33±0.05 


12.82±0.03 


11.44±0.02 


10.18±0.06 


9.63±0.06 


9.11±0.06 


8.35±0.09 


6.74±0.72 


GY92 154 


16:26:54.79 


-24:27:02.1 


M6 


20.1 


8 


17.92±0.05 


14.91±0.05 


12.87±0.05 


11.13+0.06 


10.51 ±0.06 


10.00±0.06 


9.34±0.09 


5.85±0.63 


WL21 


16:26:57.33 


-24:35:38.7 


M6 


23.8 


8 


19.0()±0.05 


15. 07 ±0.05 


12.82±0.05 


1 1.41+0. 10 


10.41 ±0.06 


9.96±0.07 


9.44±0.09 


6.93±0.65 


GY92 171 


16:26:58.41 


-24:21:30.0 


M6 


6.6 


S 


16.01±0.07 


13.11+0.02 


1 1.46±0.02 


1().09±0.06 


9.55±0.06 


9.11±0.06 


8.50±0.09 


5.80±0.63 


GY92 204 


16:27:06.60 


-24:41:48.8 


M6 


0.5 


5 


12.43±0.02 


11.40±0.02 


10.77±0.02 


9.95±0.05 


9.67±0.06 


9.32±0.07 


8.65±0.09 


5.23±0.63 


ISO-Oph 160 


16:27:37.42 


-24:17:54.9 


M6 


6.0 


5 


14.15±0.03 


12.76±0.03 


11.95±0.03 


11.00±0.05 


10.61±0.06 


10.13±0.06 


9.35±0.09 


6.74±0.63 


GY92 344 


16:27:45.78 


-24:44:53.6 


M6 


16.2 


8 


17.74±0.05 


14.8()±0.05 


12.61±0.05 


9.70±0.06 


8.71±0.06 


7.97±0.06 


7.27±0.09 


3.88±(f&3 


GY92 350 


16:27:46.29 


-24:31:41.2 


M6 


7.0 


5 


13.83±0.02 


12.21±0.02 


11.32±0.02 


10.51±0.05 


10.15±0.06 


9.73±0.06 


8.84±0.09 


6.16±(fe63 


GY92 371AB 


16:27:49.77 


-24:25:22.2 


M6 


5.4 


8 


12.78±0.02 


11.12±0.02 


10.16±0.02 


9.18±0.06 


8.76±0.06 


8.39±0.07 


7.84±0.09 


5.86±Qjj63 


GY92 397 


16:27:55.24 


-24:28:39.7 


M6 


5.0 


8 


13.04±0.02 


11.59±0.02 


10.79±0.02 


9.98±0.05 


9.62±0.06 


9.17±0.06 


8.35±0.09 


5.69±(g53 


GY92 450 


16:28:03.56 


-24:34:38.6 


M6 


20.5 


8 


19.17±0.05 


15.19±0.05 


13. 11 ±0.05 


11.72±0.05 


11.41±0.06 


11.12±0.07 


11.05±0.09 


■"A 


ISO-Oph 193 


16:28:12.72 


-24:11:35.6 


M6 


7.5 


5 


13.61±0.02 


12.02±0.02 


11.09±0.02 


10.18±0.06 


9.91±0.06 


9.43±0.07 


8.63±0.09 


5.79±0*3 














Brown dwarfs 














CD 


CFHTWIR-Oph 107 


16:28:48.71 


-24:26:31.8 


M6.25 


2.3 


this work 


14.31±0.05 


13.6()±0.05 


13.13±0.05 


12.55±0.06 


12.26±0.06 


11.88±0.07 


11.29±0.09 


9.73±(g68 


CFHTWIR-Oph 4 


16:25:32.41 


-24:34:05.2 


M6.5 


2.5 


7 


14.88±0.05 


13.93±0.05 


13.34±0.05 


12.71±0.06 


12.53±0.06 


12.37±0.11 


12.31±0.46 


... n 


CFHTWIR-Oph 106 


16:28:29.93 


-24:54:06.4 


M6.5 


4.9 


7 


15.36±0.05 


14.51 ±0.05 


13.85±0.05 


13.05±0.06 


12.57±0.06 


12.11±0.07 


11.47±0.10 


... EL 


CRBR 2322.3-1143 


16:26:23.78 


-24:18:31.4 


M6.7 


8.6 


4 


16.28±0.05 


14.62±0.05 


13.57±0.05 


12.53±0.06 


12.25±0.06 


12.02±0.08 


11.49+0.10 


8.17±f£63 


CFHTWIR-Oph 16 


16:26:18.58 


-24:29:51.6 


M6.75 


18.8 


this work 


17.09±0.05 


14.92±0.05 


13.55±0.05 


12.14±0.06 


11.61+0.06 


11.00±0.08 


10.18±0.16 


6.99±^g63 


GY92 202 


16:27:05.98 


-24:28:36.3 


M7(M6.5) 


13.0 


2,3 


16.79±0.05 


14.58±0.05 


13.15+0.05 


11.78±0.06 


11.35±0.06 


10.84±0.08 


10.19±0.09 


7.08±(g63 


CFHTWIR-OphlOl 


16:27:47.25 


-24:46:45.9 


M7 


25.8 


this work 


19.26±0.05 


16.48±0.05 


14.78±0.05 


13.35±0.06 


13.01±0.06 


12.72±0.08 


12.65±0.13 


...o 


CFHTWIR-Oph57 


16:27:04.02 


-24:02:46.9 


M7.25 


6.1 


7 


15.06±0.O5 


13.92±0.05 


13.18+0.05 


12.53±0.07 


12.37±0.07 


12.14±0.08 


11.96+0.11 


n 


CRBR 2317.3-1925 


16:26:18.82 


-24:26:10.5 


M7.5(M5.5,M7) 


10.0 


2,3,5 


14.84±0.04 


13.20±0.03 


12.14±0.02 


11.07±0.06 


10.48±0.06 


9.79±0.07 


8.74±0.09 


6.73±(J264 


CFHTWIR-Oph47 


16:26:52.26 


-24:01:46.8 


M7.5 


5.6 


7 


15.94±0.05 


14.75±0.05 


13.97±0.05 


13.35±0.07 


13.17±0.07 


13.05+0.11 


12.96±0.22 


...o 


CFHTWIR-Oph66 


16:27:14.34 


-24:31:32.0 


M7.75 


15.1 


this work 


18.42±0.05 


16.53±0.05 


15.30±0.05 


14.07±0.06 


13.49±0.06 


13.06±0.1() 


11.90=t0.11 


7.2 ±03)3 


CFHTWIR-Oph78 


16:27:26.23 


-24:19:23.1 


M7.75 


16.4 


this work 


16.35±0.05 


14.36±0.05 


13.04±0.05 


11.64±0.06 


11.14±0.06 


10.62±0.07 


9.93±0.09 


7.79±Qa63 


CFHTWIR-Oph96 


16:27:40.84 


-24:29:00.8 


M7.75 


2.4 


this work 


14.60±O.O5 


13.76±0.05 


13.19+0.05 


12.51±0.06 


12.24±0.06 


11.79±0.07 


10.96±0.09 


8.33±(to64 


GY92 3 


16:26:21.90 


-24:44:39.8 


M8(M7.5) 


0.0 


6,5 


12.34±0.02 


11.48±0.02 


10.86±0.02 


9.97 ±0.06 


9.57±0.06 


9.21±0.06 


8.46±0.09 


5.82±(|63 


GY92 64 


16:26:32.53 


-24:26:35.4 


M8 


11.0 


2 


16.36±0.05 


14.44±0.05 


13.17±0.05 


12.15±0.06 


11.94±0.06 


11.71+0.08 


11.77±0.17 


... 3 


GY92 264 


16:27:26.58 


-24:25:54.4 


M8 


0.0 


6 


13.00±0.02 


12.35±0.02 


11.84±0.02 


11.09±0.06 


10.61±0.06 


10.18±0.07 


9.32±0.09 


6.36±0|63 


CFHTWIR-Oph34 


16:26:39.92 


-24:22:33.6 


M8.25 


9.7 


7 


15.97 ±0.05 


14.53±0.05 


13.48±0.05 


12.21±0.07 


11.66±0.06 


11.13±0.24 


10.16±0.40 


... £ 


GY92 11 


16:26:22.27 


-24:24:07.1 


M8.5(M6.5,M8.5) 


8.0 


9,2,5 


16.18±0.05 


14.73±0.05 


13.75±0.05 


12.39±0.06 


11.56±0.06 


10.79±0.07 


9.61±0.09 


6.12+0^54 

... CD 


GY92 141 


16:26:51.28 


-24:32:42.0 


M8.5(M8) 


0.7 


1,4 


15.06±0.05 


14.35±0.05 


13.86±0.05 


13.15+0.05 


12.84±0.06 


12.33±0.07 


11.78+0.10 


GY92 310 


16:27:38.63 


-24:38:39.2 


M8.5(M7,M6) 


5.7 


2,3,5 


13.27±0.02 


11.93±0.02 


11.08±0.02 


10.18±0.06 


9.77±0.06 


9.34±0.06 


8.35±0.09 


4.78±(§63 


SONYC-RhoOph-1 


16:26:56.33 


-24:42:37.8 


M9 


5.0 


10 


17.46±0.05 


16.44±0.05 


15.73±0.05 


14.74±0.06 


14.46±0.07 


14.22±0.23 






CFHTWIR-Oph77 


16:27:25.64 


-24:37:28.6 


M9.75 


10.0 


this work 


18.23±0.05 


16.54±0.05 


15.36±0.05 


14.16±0.06 


13.9 ±0.06 


13.52±0.12 


13.42±0.20 


'"B" 


CFHTWIR-Oph98 


16:27:44.20 


-23:58:52.4 


M9.75 


3.0 


this work 


17.08±0.05 


15.78±0.05 


14.98±0.05 


13.66±0.08 


13.32±0.07 


12.98+0.11 


12.28+0.13 




CFHTWIR-Oph9 


16:26:03.28 


-24:30:25.9 


LO 


6.6 


this work 


17.76±0.05 


16.33±0.05 


15.32±0.05 


14.12±0.07 


13.63±0.07 






...o 


CFHTWIR-Oph 18 


16:26:19.41 


-24:27:43.9 


LO 


7.8 


this work 


18.92±0.05 


17.16±0.05 


16.04±0.05 


14.68±0.03 






11.73±0.16 


■ ■■3- 


CFHTWIR-Oph90 


16:27:36.59 


-24:51:36.1 


LO 


2.4 


this work 


16.83±0.05 


15.65±0.05 


14.85±0.05 


13.86±0.06 


13.53±0.06 


13.15±0.08 


12.32±0.09 


8.91+(E68 


CFHTWIR-OphlOO 


16:27:46.54 


-24:05:59.2 


LO 


8.0 


this work 


17.94±0.05 


16.34±0.05 


15.26±0.05 


14.41 ±0.07 


14.17±0.08 


13.96+0.15 




o 

■■■3'. 


CFHTWIR-Oph 103 


16:28:10.46 


-24:24:20.4 


L0 


7.5 


this work 


17.74±0.05 


16.15±0.05 


15.07±0.05 


13.98±0.06 


13.68±0.06 


13.29+0.11 


12.97±0.14 




CFHTWIR-Oph33 


16:26:39.69 


-24:22:06.2 


L4 


3.6 


this work 


18.16±0.05 


16.74±0.05 


15.68±0.05 


14.35±0.10 


14.12±0.08 









{a) Spectral types determined from the following st udies: l.lLuhman et alj<1997l) : 2. IWilking et all J1999I): 3. lLuhman & RiekeUl9991) : 4. ICushing et"aI1 <2000b ; 5. lNatta et all j2002h : 6. IWilking et"aT]<2005l) : 

7. lAlves de Oliveira et aljd2010h : 8. lMcClure et all 120101) : 9. IComeron et alJfcOld) : lO. lGeers et alJfeOllh . 



C. Alves de Oliveira et al,: Spectroscopy of new brown dwarf members in p Ophiuchi 



Table[5]lists the 29 brown dwarfs and 13 very low-mass stars 
at the substellar limit members of the cluster with spectroscopic 
confirmation known to date with coordinates and near-IR pho- 
tometry, which we adopt in the following sections to study their 
properties. The JHK S photometry prov ided for each member 
is fro m the WIRCam/CFHT survey dAlves de Oliveira et al.l 
1201 Ol) . except for members that were saturated or out of the 
field of view, for which we have listed 2MASS photometry. For 
each member, the spectral type and Ay are shown as compiled 
from the literature or derived from the numerical spectral fitting 
in this work. We have mid-IR photometry from t he Spitzer 
space telescope, taken by the C2D legacy project (lEvans et al.l 
I2003L From Cores to Disks), which mapped the p Ophiuchi 
molec ular cloud with Spitzer' 's Infrared Camera (Fazi o et alj 
12004 IRAQ in the 3.6, 4.5, 5.8, and 8.0 pm bands over a 
region of 8.0 deg 2 and with the Multiband Imaging Camera 
dRieke et al.l 12004 MIPS) in the 24 and 70 fan bands over a 
total of 14.0 deg 2 dPadgett et alj|200 8;). The data were retrieved 
from the C2D point - source catalogues of the final data delivery 
( lEvans & c2d Team! 120051) . using the NASA/ IPAC Infrared 
Science Archive. All fluxes were converted to magnitudes using 
the following zero points: 280.9+4.1, 179.7+2.6, 115.0+1.7, 
64.1+0.94 (Jy), for the 3.6, 4.5, 5.8, and 8.0 fim IRAC bands, 
respectively, and 7.17+0.11 (Jy) for the 24 fan MIPS band. 
The Spitzer catalogues were merged with the WIRCam or 
2MASS catalogue, requiring the closest match to be within 1". 
A counterpart was found for all substellar members, which were 
detected in one or more mid-IR bands. Only detections above 
2 cr were kept, though all objects have 7 cr detections in the 
[3.6] and [4.5] bands. 



5.2. Slopes of the spectral energy distributions 

The disk population of the p Ophiuchi cluster has been stud- 
ied in great detail (e.g.. [Natta et al.ll2.002i: [A ndrews & Williams 
20071 iFurlan et al]l2009t lAndrews et al|2010tlRicci et alJ l2010: 
Cieza et al.l 120101: iMcClure et all I20101) . It is found that, de- 
spite the young age of the cluster (~lMyr), many circumstellar 



Table 6. SED slopes for spectroscopically confirmed brown 
dwarfs in p Oph. 



disks in p Oph are evolved showing feature s such as mi llime- 
ter size dust grains in the outer disk r egion (iFurlan et al] l"2009: 
lAndrews et alj20iotlRicci et alj201(i|). inner o pacity holes char- 
acteristic of transition disks (Cieza et al.l2010 ). dust settling, and 
dust grain processing (e.g.. IMcClure et aU 120101) . In this sec- 
tion, and in light of the newly discovered substellar members 
in our survey, we use SEDs to classify the evolutionary state 
of each object. The slope of the SED {a=dlog{AF / \)/dlog{A)) is 
used to distinguish YSOs into various cl asses according to their 
mid-IR (and millimetre) emittin g energy dLada & Wilk in g 1984; 
Ada ms et alJI 1987k lLadall 19871). re presenting deeply embedded 
YSOs (Class 0. lAndre et alj|1993b . embedded YSOs still sur- 
rounded by an envelope (Class I), YSOs possessing a circum- 
stellar disk (Class II), flat-spectrum sources which may represent 
an int ermediate stage between Class I and II (FS, iGreene et alj 
11994 . and sources with a weak or debris disk (Class III). 

Following the approach of disk studies using IRAC and 
MIPS phot ometry in other young clus ters (e.g., Chamaeleon I 
and Taurus, Luhman et aDl2008aLl2010h . we computed the spec- 
tral slopes across several bands (2sT s /[8.0], K S /[2A], [3.6]/[8.0], 
[3.6]/[24]), which allowed us to derive an SED class even for 
members which were not detected at 24 pm. All fluxes were first 
dereddened by the respective extinction value lis ted in Table |5j 
using the extinction law from Flah erty et al.l d2007h . We have de- 



ID 


Ks/[8.0] 


Ks/[24] 


[3.6]/[8.0] 


[3.6]/[24] 


Class 


GY92 15 


-1.09 


-1.07 


-0.84 


-0.95 


n 


GY92 109 


-1.20 


-1.41 


-1.04 


-1.40 


11 


GY92 154 


-1.18 


-0.68 


-1.23 


-0.57 


11 


WL21 


-1.44 


-1.20 


-1.11 


-1.00 


11 


GY92 171 


-0.99 


-0.88 


-1.16 


-0.92 


11 


GY92 204 


-1.32 


-0.78 


-1.36 


-0.65 


11 


ISO-Oph 160 


-1.22 


-1.03 


-1.08 


-0.93 


11 


GY92 344 


0.27 


0.07 


-0.41 


-0.27 


11 


GY92 350 


-1.35 


-1.08 


-1.08 


-0.89 


11 


GY92 371AB 


-1.39 


-1.37 


-1.42 


-1.38 


11 


GY92 397 


-1.29 


-1.05 


-1.08 


-0.90 


11 


GY92 450 


-2.23 




-2.54 




in 


ISO-Oph 193 


-1.38 


-1.03 


-1.23 


-0.88 


11 


CFHTWIR-Oph 107 


-1.59 


-1.64 


-1.45 


-1.59 


n 


CFHTWIR-Oph 4 


-2.17 




-2.44 




in 


CFHTWIR-Oph 106 


-1.33 




-1.14 




11 


CRBR 2322.3-1143 


-1.70 


-1.02 


-1.84 


-0.90 


11 


CFHTWIR-Oph 16 


-1.23 


-0.82 


-1.01 


-0.62 


11 


GY92 202 


-1.27 


-0.87 


-1.31 


-0.78 


11 


CFHTWIR-Oph 101 


-2.41 




-2.62 




in 


CFHTWIR-Oph 57 


-2.19 




-2.33 




in 


CRBR 2317.3-1925 


-0.83 


-1.05 


-0.39 


-0.93 


11 


CFHTWIR-Oph 47 


-2.32 




-2.53 




hi 


CFHTWIR-Oph 66 


-1.05 


-0.15 


-0.68 


0.25 


11 


CFHTWIR-Oph 78 


-1.31 


-1.27 


-1.24 


-1.23 


11 


GY92 3 


-1.10 


-0.95 


-1.11 


-0.92 


11 


GY92 64 


-2.28 




-2.66 




in 


GY92 264 


-1.02 


-0.79 


-0.81 


-0.64 


11 


CFHTWIR-Oph 34 


-0.87 




-0.70 




11 


CFHTWIR-Oph 96 


-1.27 


-1.05 


-1.09 


-0.92 


11 


GY92 11 


-0.22 


-0.16 


0.18 


0.03 


11 


GY92 141 


-1.36 




-1.29 




11 


GY92 310 


-1.11 


-0.61 


-0.87 


-0.37 


11 


SONYC-RhoOph-1 










m? 


CFHTWIR-Oph 77 


-1.86 




-2.22 




in 


CFHTWIR-Oph 98 


-1.02 




-1.33 




11 


CFHTWIR-Oph 9 










n? 


CFHTWIR-Oph 18 


-0.09 




0.38 




11? 


CFHTWIR-Oph 90 


-1.11 


-0.67 


-1.13 


-0.55 


11 


CFHTWIR-Oph 100 










m? 


CFHTWIR-Oph 103 


-1.6 




-1.85 




in 


CFHTWIR-Oph 33 













rived in this manner the spectral slopes for 38 of the low-mass 
stars at the substellar limit and brown dwarf members. 



5.3. Disk evolutionary states 

Based on the spectral slopes derived, we classify each object into 
Class II or Class III. Figure [6] shows the histograms of the slope 
values (upper panel) and the distribution of slopes per spectral 
type (lower panel), which is compared to the slo pes of the pho- 
tosph eric colours for young objects (filled line, Luhm an et al.l 
2010). For a(KJ[8.0]) and a([3.6]/[8.0]) the histogram shows 
two well-separated groups corresponding to Class II (> -2) and 
Class III (< -2) objects. The same groups are also clearly sep- 
arated in the plot of the slopes as a function of the spectral 
type, where all objects classified as Class III follow the line 
representing stellar photospheres, while the Class II have sys- 
tematically higher values. The only exception is CFHTWIR- 
Oph 103, classified as a Class III since it shows a slope that 
is only marginally higher then expected for the photospheric 
colours of a young L0. Due to the limited sensitivity of the 24 pm 
measurements, none of the Class III objects has a detection at 
this wavelength, and therefore does not appear in the a(K s /[24]) 
and ff([3.6]/[24]) plots. These diagrams recover, however, all the 
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M6 M7 M8 M9 LO 
Spectral Type 



M6 M7 M8 M9 LO 
Spectral Type 



M6 M7 M8 M9 LO 
Spectral Type 



M6 M7 M8 M9 LO 
Spectral Type 



Fig. 6. SED slopes for spectroscopically confirmed brown dwarfs in p Oph. The histograms of the different spectral slopes are shown 
in the upper panel. SED slopes versus s pectral type are dis play in the lower panel, where the solid line represents the slopes of the 
photospheric colours for young objects (Luhman et al. 2010h. 



Class II. Inside the Class II group are two objects that have SEDs 
indicative of Class I ( g ~0), GY92 344 and CFHTWIR-Oph 66. 
dMcClure et al.l2010h find evidence for a disk around GY92 344, 
so we adopt this result and classified it as Class II. CFHTWIR- 
Oph 66 shows only moderate excess from 3.6 to 8.0 pm, and is 
more consistent with a Class II. 

Four objects that are not detected at 8 pm or 24 pm (SONYC- 
RhoOph-1, CFHTWIR-Oph 9, 100, 33) nevertheless have de- 
tections at other IRAC bands, which we used to constraint their 
SED class. CFHTWIR-Oph 100 and SONYC-RhoOph-1 are de- 
tected in the three other IRAC channels. In the IRAC colour- 
colour diagram [3.6]-[4.5] vs. [4.5]-[5.8], both objects occupy 
the same region of the Class III objects, not showing a signif- 
icant excess in any of the colours. Based on this evidence, we 
classified them as likely Class III (denoted as III?). CFHTWIR- 
Oph 9 and 33 are only detected at 3.6 and 4.5 pm. In the colour- 
colour diagram K s -[3. 6] vs. [3.6]-[4.5], CFHTWIR-Oph 9 shows 
a clear excess mid-IR excess, so we classify it as likely Class II 
(denoted as II?). CFHTWIR-Oph 33 occupies a region between 
Classes II and III, so we cannot ascertain its SED class. Finally, 
CFHTWIR-Oph 18 is detected at 3.6 pm but has a 2<x detection 
at 8.0 pm. The slope of the SED shows it has strong mid-IR ex- 
cess, which places it in the Class I/Class II range, though given 
the limited information at mid-IR wavelengths this classification 
remains provisory (denoted as II?). 

According to the SED classifications, 68+j % of the brown 
dwarf members of p Oph are Class II sources, i.e., have an 
SED representative of a disk. Statisti cal errors are calcu lated 
following the method proposed in Burga sser et al.l (120031) . If 
we also include all objects at the substellar limit (spectral type 
~M6), the percentage of Class II sources goes up to 76*5%. The 
disk frequency in the stellar domain of p Oph has been esti- 



mated as 27+5% for an extinction-limited sampl e of Ay<8 and 
repre sentative of YSOs with masses >0.2 Mq (Erickson et al. 
2011). We considered the substellar members that lay within 
the same extinction limit (see Fig. |9j, and derived a disk frac- 
tion of 74%^j 2 %, confirm ing the large number of brown dwarfs 
with disks in this cluster. iLuhman et al.l (1201 Oh have found the 
disk fraction in Taurus (~1 Myr) to be 75% for objects with 
>0.3 Mq, and 45% for low-mass stars and brown dwarfs. By 
comparing these results to the disk fraction in IC 348, Cha I, 
and o~ Ori, the authors suggest stars with masses greater than 
0.5 Mq members of dense clusters have shorter disk lifetimes. 
Th e results of the disk fr action of the stellar population derived 
bv lErickson et al.1 (1201 ll) also support this hypothesis. The disk 
fraction we derive for brown dwarfs (> M6) is higher than wha t 
is found, for example, in cr Ori (60%. ILuhman et alJ [2008b), 
IC 34 8 (51%.lLada et al.ll2006llLuhman et al.ll2008ah . and Cha I 
(44%. ILuhman et al.l2008al) . A detailed analysis of the disk frac- 
tion per spectral type for the p Ophiuchi cluster across the entire 
mass range is needed to understand this result, though such a 
study is beyond the scope of this paper. 



6. The Rho Ophiuchi population 

To compare the properties of the population of the p Ophiuchi 
molecular cloud to other young clusters, it is fundamental to 
analyse this newly found population of substellar members in 
respect to the stellar content of the cluster. Despite being one 
of the closest star-forming regions, the high and variable extinc- 
tion in p Ophiuchi and the fact that the central cluster extends 
over an area on sky greater than 1 deg 2 means that many can- 
didate members still lack a spectroscopic classification, hence a 
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reliable membership assessment. Here we summarise the infor- 
mation we gathered from different spectroscopic surveys. 

6.1. The cluster's stellar population 

We took as a starting point for the membersh i p com pilation the 
review done on the cluster by Wilking et al. (2008), excluding 
the substellar members and members with a spectral type of M6, 
which are at the substellar limit and have already been discussed 
in the previous section. To this list, we added newly found mem- 
bers with spectroscopic confirmation published afterwards. The 
literature review can be summarised as follows: 



- |Wilking et al.l(l2008l) : in total, the compiled list contains 316 
objects, from which 27 are already included in Table[5] From 
the remaining 289 listed targets, 156 members have assigned 
spectral types from optical and/or near-IR spectroscopy and 
evidence of youth (please see the review for references), 
while 133 did not have any spectral type assigned at that 
time; 

- ICieza et ail (120101) derived spectral types for a total of 34 
sources, nine of these had previously determined spectral 
types, while the remaining 25 sources are located outside the 
central region of the cluster and were therefore not included 
in our compilation. 

- lAlves de Oliveira et alJ (|2010) computed spectral types for 
six new brown dwarf s and one very low mass star. 

- iMcClure et al.l (|2010) determined spectral types for a few 
tenths of sources and their binary components (from which 9 
are show n in Tabled, includin g 25 candidate members com- 
piled by Wilking et al. (2008) which lacked spectroscopic 
and four sources that were not part of that list but have spec- 
tral types and membership determined in the literature. N.B., 
some of the s e spec tral types had been published earlier by 
iFurlan et alJ ([2009), though it is in this paper that the data 
reduction and analy sis are presented). 

- iGeers et all d201 lb derived the spectral type for one new 
brown dwarf. 

- lErickson et alJ (1201 ll) confirmed 30 new members with op- 
tical spectral types, where only three were part of the 
IWilking et alJ d2008h compilation. They combined all the 
new members with 88 previously known members with 
optical spectral typ e from an analogous optical survey 
(IWilking et al.ll2005h . 14 members compiled from the liter- 
ature, and three members for which they presented optical 
spectral t ypes but had already nea r-IR spectral types deter- 
mined bv lLuhman & Riekd dl999h . 

- In this work we confirm spectroscopically three new mem- 
bers with spectral types indicative of masses above the sub- 
stellar limit (besides the 13 brown dwarfs already included in 
Table |5j and the four members with uncertain spectral types 
that we do not include in this analysis); 

totalling a list of 250 members of the cluster where 208 have 
spectral types earlier than M6, and 42 have spectral types later 
than or equal to M6. 

Several other so urces have been associated with the 
cloud through X-rays dlmanishi et al] 1200 It : iGagne et al.| [2004; 
lOzawa et all 12005b [PilUtteri et alJ l2010h . near-I R variability 
dAlves de Oliveira et al.1 1201 Oh. mid-IR excess (Padgett ~et al] 



magnitudes, e xtrag alactic contamination is present. In the 
Wilk ing et al] d2008l) list of candidate members, there are 51 
sources that have been associated with the cloud from X-ray sur- 
veys but not confirmed spectroscopically. From these, only 1 1 
candidates are within our photometric completeness limits, and 
more conservative extinction limited sample (see Sect. 16.2. j} , 
with //-band magnitudes ranging from 9.5 to 14 mag. We there- 
fore do not expect our results to be affected in a significant 
way. Candidate samples found by the other aforementioned tech- 
niques are known to have much higher contamination rates and 
are therefore not included in the following sections. We note, 
however, that the majority of these candidates are outside our 
most conservative extinction limited sample, as well as our com- 
pleteness limits (see Sect. [6T27T]). 

While this paper was in the process of revision. iMuzic et"ail 

d201 ll) presented spectra for 21 members of p Oph, 14 of which 
had spectral types already published in the literature and five 
are simultaneously discovered in this paper (CFHTWIR-Oph 9, 
16, 31, 58, 78). The only objects not included in our census are 
CFHTWIR-Oph 58 for which we could not assign a spectral type 
(corresponds to SONYC-RhoOph- 6 with a spectral typ e of M6.9 
calculated from the H2O defined in lAllers et al.l ([2007)), and two 
low-mass stars (SONYC-RhoOph-3 and 4), which have visual 
extinctions greater than 8 ma g (see Se ct. 16.2.11 ). Therefore, we 
conclude that the results from Muzicet aDd2011) do not affect 
any of the conclusions of this paper. 



6.2. The H-R diagram 

To study the masses of the members of the cluster, we first place 
all the members spectroscopically confirmed in the Hertzsprung- 
Russell (H-R) diagram (Fig. [7]). For the 208 stellar members 
of the cluster, we compiled information from the literature on 
the spectral type, Ay, and near-IR magnitudes. Fourteen ob- 
jects did not have any measured extinction and are not be dis- 
played in the H-R diagram but are included in the other sec- 
tions. To convert spectral types to temperature, we adopted the 
temperature scale from Sch midt-Kahler (1982) for stars earlier 
than M0, and the scale from iLuhman et al] d2003bl) for sources 
with spectral type between M0 a nd M9.5. For the L dwarfs, 
we applied the scale proposed by lLodieu et al] d2008l) extrap- 
olated to the L4 spectral type. The temperatures derived in this 
way have an estimated uncertainty of +150 K. The bolometric 
luminosity was calculated using the dereddened J magnitude, 
a distance to the cloud of 13 pc, and bolometric correc tions 
from Schm i dt-Ka hler (1982), lKenvon & Hartmairnl d 1995b . and 
iDahn et all d2002l) . For the most massive members with spectral 
type B, we plotted the lumino sities derived from optical photom- 
etry by Ericks on et al] d201 lb . Figure [7] shows all the members 



overplotted with NextGen and Dusty isochrones (Baraffe et al 



2008; Wilking et al] 120081; iGutermuth et al l 120091; iMarsh et al] 
2010), and methane photometric colours (Hais ch et al.l 12010). 
but lack spectroscopic confirmation. X-ray surveys have low 
contamination rates, particularly when large amounts of ex- 
tinction block background contaminants, though at faint 



1998; Chabrie retani2000l) and the isochrones from Siessetal 
(120001) . for 1, 3, 10, and 30 Myr. 

As seen in other young clusters and associations, a large ap- 
parent spread in ages is seen in the H-R diagram, with most of 
the members lying between the 1 Myr and 10 Myrs isochrones, 
and some objects already close to the 30 Myr i sochrone. The age 
for p Oph has been estimated as ~0.3-l Myr dGreene & Meverl 
ll995HLuhman & Riekdfl999l) bv comparing the location on the 
H-R diagram of near-IR spectroscopically confirmed members 
found in the central core of the L1688 cloud with pre-main- 
sequence evolutionary models. The proximity on the sky be- 
tween the p Oph cluster and the southeastern edge of the Upper 
Sco association (which adjoins the Sco OB2 association), has 
been several times used to explain the high star-forming activ- 
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Fig. 7. H-R diagram for all the spectroscopically confirmed 
members of p O phiuchi. The models are the NextGen and 
Dusty isochrones dBaraffe et al1ll998t lChabrier et al.ll2000h (up- 
per panel) and the isochrones from ISiess et aL ( 20001) (lower 
panel) for 1, 5, 10, and 30 Myr, labelled with mass in units of 
M . 



ity seen in this cluster. That a supernova explosion occurred 
~ 1 .5 Myr ago in Upper Sco is usually invoked to explain a sce- 
nario where the strong s hock wave "triggered" star formatio n 
in the nearby p Ophiuchi (lLoren & Woottenl[T986tlLorenl[l 989). 
matching the ages derived for its YSO population. However, this 
somehow simplistic interpretation should be reviewed in light of 
more recent theories of cloud formation and lifetime (see, for 
exam ple, the discussion on Sco-Cen by iPreibisch & Mam aiek 
2008). Based on several observations, one can nevertheless ex- 
pect that the on-going star formation in p Oph is i nfluenced by 
the ne arby high-mass members of Upper Sco (e.g.. lNutter et al.l 
2006). Observations of a low extinction population, also com- 
monly called the surface population, used optical spectroscopy 
to confirm many members and by placing then on the H-R 
diagram derive an older age of 1 to 5 Myr for these objects 
( Bou vier & Appenzellerll 19921: iMartin et al.lll998r.lWilking et al l 
120051) . This has been interpreted bv lWilking et al.l l 2005 ) as an 
indication that star-formation might have started even earlier in 
Ophiuchus, simultaneously with Upper Sco. However, there is 
an increasing scrutiny of the interpretation of the luminosity dis- 
persion seen i n the H-R diagram as a true spread in the age of 
the stars (e.g., Jeffr ies et al.ll201 lh . and therefore it could be er- 
roneous to formulate a star formation scenario based only on 
the luminosity spread. For example, several morphological prop- 
erties of young stellar objects can affect its derived luminosity, 



such as the presence of circumstellar emission, variability caused 
by photospheric spots or the inner disk, and its accretion history. 
Observationally, unresolved binary systems will also result in an 
incorrect luminosity, objects observed edge-on, as well as errors 
in determining extinction, which is non-negligible in this cluster. 
Finally, errors in the spectral type classification and the fact that 
we adopt the same distance to all stars will also significantly af- 
fect the position in the H-R diagram. In Fig. [7] the few objects 
located below the 30 Myr evolutionary track have large variabil- 
ity (Oph 30, see Sect. I4.1.2I I. while the others are Class I and 
Class II sources, where the circumstellar material is most likely 
compromising one or more parameters used in deriving their lu- 
minosity. 



6.2.1 . Spatially and extinction limited sample 

To ensure that the properties we derive for the p Ophiuchi clus- 
ter are representative of its population, a sample that is spa- 
tially and extinction limited, while maximizing completeness 
must be constructed. In a first step, we select only targets that 
are contained in the WIRCam survey field (Fig. [8), since we 
have explored the substellar population of the cluster over this 
area in an unbiased way. All the selected members are plotted 
on the colour-magnitude diagram H vs. J-H (Fig. [9), where the 
comple teness limit shown is that of the WIRCam photometric 
survey (Al ves de Oliveira et a n i2oloh . In this diagram, we also 
plot the 1 Myr isochrone at a distance of 1 30 pc from the BT- 
Settl evolutionary model (lAllard et alJl20ld) . which ensures a 
smooth transition between the stellar and substellar reg imes. The 
1 Myr isochrones for the NextGen dBa raffe et ai] |1998l) . Dusty 
(Chab rieretaTl2 000). and lSiess et al .1 (120001) are also shown for 
comparison. The isochrone has been reddened by 8, 15, and 
20 mag of Ay, encompassing different depths of the cloud. The 
extinction-limited sample with Ay<8 is likely to be complete at 
all mass ranges, since the combination of optical and near-IR 
spectroscopic surveys in the stellar regime have targeted these 
objects. In the substellar regime, we observed all but one candi- 
date brown dwarf (CFHTWIR-Oph 72) spectroscopically down 
to the completeness limit, defining in this way a nearly com- 
plete sample. The limits of Ay of 15 and 20 mag are probably 
not complete. In the substellar regime, the Ay<15 limit contains 
two brown dwarf candidates that have not been observed spec- 
troscopically, and this number rises to three if we go to <20. 
However, the incompleteness is greater in the stellar regime, 
where at this higher extinction regions, several candidates from 
X-ray surveys still lack a spectroscopic confirmation and have 
therefore not been included (see Sect. 16. U . Hereafter, we study 
the properties of the population with Ay<8, which we consider 
to be representative of the cluster, but also display the results for 
higher extinction regions with the caution that these likely suffer 
from incompleteness in the stellar regime. 



6.3. The IMF 

In this section, we analyse the distribution of luminosities and 
masses of the p Oph population and compare it to other star- 
forming regions. Figure \W\ shows the spectral type distribution 
for our extinction-limited sample. The peak occurs for a spec- 
tral type M5 (bin contains members with spectral types >4.5 
and <5.5), which is similar to the spectral type distributions 
of other star-formin g regions such as IC 348 and Cha I (e.g., 
iLuhman et al.l 12009). The distribution of the deredden H mag- 
nitudes corrected for a distance of 130 pc is shown in Fig. Q~T] 
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Fig. 8. Spatial distribution of all the spectroscopically confirmed 
members of p Ophiuchi as compiled from the literature and de- 
termined in this paper. The diagram shows the stellar members 
(black filled circles) and members at the substellar limit and 
brown dwarfs (grey filled circles) , superposed on th e contours of 
the COMPLETE extinction map ( (Ridge et al.l2006l) . The field of 
the WIRCam survey is shown and is used to define a spatially- 
limited sample. 



shows a broad peak at an absolute magnitude of 4 to 5. We con- 
structed the IMF for the spectroscopically confirmed members of 
p Oph by deriving masses from the 1 Myr evolutionary models, 
according to each target's effective temperature. We did not in- 
clude members with uncertain spectral types, and close binaries 
that are not resolved in 2MASS or WIRCam images are treated 
as a single object. The resulting IMF is shown in Fig. [12] and it 
exhibits a peak at ~0.12 M Q . This result also agrees with what 
is found for other young clusters and the field mass function, 
which have ma ss function peakin g in the 0. 1 -0.2 M G range (see, 
for example, Bastian et al. 2010, for a review). 

More quantitatively, we can calculate the number of brown 
dwarf to stars, defined as R = N(0.02 < M/M Q < 
0.08)/JV(0.08 < M/M < 10). This ratio is of 0.09+J}^ for 
the A v <8 sample (and 0.09+° °2 for a n mem bers, 0.13+°-^ for 
A v <15, and 0.13+°jg for Ay<20). This value is lower than 
what has been found in other 1-2 Myr old clusters (0.13 in 
IC 348, 0.18 in Taurus, 0.26 in Cha I and 44 in NGC 1333 
iLuhman et al.ll2003br. lLuhmadl2004b[ 120071: IScholz et al.ll20TTl 
respectively). However, this ratio is highly dependent on the 
accuracy of the spectral classification at the substellar limit. 
As pointed out in Sect l5.11 13 members of the cluster have 
been classified with an M6 spectral type through near-IR spec- 
troscopy, a technique that usually carries an error +1 subclass. 
According to the evolutionary models used to determine the 
masses, all these objects have masses greater than ~0.08 M . A 
spectral type of ~M6.25 would, however, place these objects in 
the substellar regime and significantly change the ratio of brown 
dwarfs to stars. The value derived for this ratio, if we place the 
members classified as M6 in the substellar domain, would then 



5+0.04 



be 0.14+;j{£ for A v <8 (and 0.16+;j ^ for all members, 0.19! 0Q , 
for A y <15, and 0.20+°^ for A v <20). These values agree with 



Fig. 9. Colour-magnitude diagram for the spatially-limited sam- 
ple of members of p Ophiuchi with a known spectral type. The 
1 Myr isochrone shif ted at a distance of 1 30 pc from the BT-Settl 
evolutionary model (lAllard et al.ll20ldl) is shown for Ay=0, 8, 
15, and 20 mag (solid black line), while th e 1 Myr isochrones fo r 
the NextGen fearaffe et aljfl9"98l) . Dusty dChabrier et all llzOOO), 
and Siess et al .1 (120001) are also shown for comparison (solid grey 
line). The diagram shows the stellar members (black filled cir- 
cles), members at the substellar limit and brown dwarfs (grey 
filled circles), members with uncertain spectral types (open cir- 
cles), contaminants found in our spectroscopic follow-up (as- 
terisks), and candidate brown dwarfs from the WIRCam survey 
which have not been observed spectroscopically (crosses). The 
dashed line shows the photometric completeness limits of the 
WIRCam survey at 7=20.5 and //=18.9. 
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Fig. 10. Histogram of the spectral types for the population of 
the p Ophiuchi cluster. The lines show the distribution of all 
members (dotted line), and for the extinction limited samples 
for Ay <8 (solid line), 15 (dashed line), and 20 mag (dash-dotted 
line). The peak of the distribution is located at a spectral type of 
M5. 
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Fig. 11. Histogram of the absolute H magnitudes corrected for 
extinction and assuming a distance of 130 pc. The distributions 
are depicted in the same linestyles as Fig. [TO] The histogram 
peaks at a magnitude of 4. 
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Fig. 12. IMF for the spectroscopically confirmed population of 
the p Ophiuchi cluster. The peak of the distribution is located at 
-0.12 Mq. Linestyles as in Fig.fTOl 



those found for the aforementioned young clusters, with the ex- 
ceptionofNGC 133 3 claimed to have an excess of brown dwarfs 
dScholz et al.ll201 il) . Given that the errors associated with the 
spectral type classification imply an error in the ratio of brown 
dwarf to stars which is higher than the purely statistically derived 
errors, we conclude that within these limitations, our findings 
agree with other young clusters, and support neither the scenari o 
of a paucity of brown dwarfs in the cluster (Erick son et al.l20TTb . 
nor an excess of brown dwarfs down to 0.02 M . 

7. Conclusions 

We have conducted a spectroscopic follow-up of candidate 
members of p Ophiuchi uncovered in our photometric survey. 
We combine the results of this study with a census of the clus- 
ter, which was compiled from the literature to provide an up- 
dated study of the YSO population of the cluster. Our findings 
are summarised as follows. 



- The spectroscopic follow-up resulted in the confirmation of 
20 new members of the cluster, from which 13 are brown 
dwarfs, three are low-mass stars, and four have uncertain 
spectral types. Combined with our previous spectroscopic 
follow-up, this amounts to 19 new brown dwarfs members 
of p Oph found through the WIRCam survey. Fourteen of 
the observed candidates are classified as likely contaminants. 
Within the photometric completeness limits, we find the suc- 
cess rate of the WIRCam survey to be between 49% and 
74%. 

- From the new brown dwarfs uncovered in this study, four 
have a spectral type of L0 and one brown dwarf is an L4. 
These are the first L-type objects to be found in the clus- 
ter, and according to evolutionary models, these objects have 
masses between ~4 and ~10 Jupiter masses, placing them 
amongst the least massive young brown dwarfs known to 
date. 

- We compiled a list of confirmed brown dwarfs in the cluster 
and using Spitzer mid-IR photometry, we derived a ratio of 
disks of 76^%. This value is higher than what is found in 
other young clusters, so further studies are needed to under- 
stand the implications of these result. 

- A census of the spectroscopically confirmed population of 
p Ophiuchi was compiled from the literature, totalling 250 
members. By constructing a spatially and extinction-limited 
sample, which we believe is representative of the cluster's 
population, we looked at the luminosity and mass distribu- 
tions. We find that the distribution of spectral types peaks at 
M5, similar to what is seen in IC 348 and Cha I. The mass 
function peaks at 0.12 Mq, which is consistent with what is 
found in other clusters and the field. Taking the errors associ- 
ated with the spectral type classification into account, we find 
a ratio of brown dwarfs to stars that also in agrees with the 
one derived for other young clusters. We therefore conclude 
that there is no evidence for variations in the mass function 
of p Oph when compared to other nearby young clusters. 
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